Comment on "Absence of a Slater Transition in 
the Two-Dimensional Hubbard Model" 

While we agree with the numerical results of Ref . [1] , 
we arrive at different conclusions: The apparent open- 
ing of a gap at finite-temperature in the two-dimensional 
weak- coupling Hubbard model at half-filling. does not ne- 
cessitate an infinite correlation length £ (Slater mecha- 
nism) nor a thermodynamic finite-temperature metal in- 
sulator transition (MIT). The pseudogap is a crossover 
phenomenon due to critical fluctuations in two dimen- 
sions, namely to the effect of a (tt, 7r) spin-density-wave 
(SDW) £ that is large compared with the thermal length. 

We use the units of Ref. [1]. The inset of Fig.l 
shows (n^ni) obtained in Ref. [1] for N c = 36, 
U = 1 and N c = 64, U = 0.5 along with 
the corresponding results obtained [2] from the lo- 
cal moment sum rule (T/N c ) J2 q Xsp (?) = 1 — 
2 (n^ni), supplemented with the relations xlp 1 (?) — 

Xoiqy 1 - % and U sp = U (n T 7ij> / ((ra T > (n^». 
Fig.l also shows the pseudogap in the density of 
states p(u>) obtained from [3] (&) = f/n_ CT + 
TTF E 9 i3U sp Xsp(q) + U chXch (q)\ Gl(k + q) which in- 
eludes the effects of both spin Xsp and charge Xch 
fluctuations and satisfies ^Trp^G"] = J7 (n^n;) .The 
charge fluctuations are constrained by the sum rule 
(T/N c ) J2 q (Xsp (?) + Xch (?)) = 1- As temperature is 
lowered from T = 1/20 to 1/22 and 1/32, the pseu- 
dogap in p (u>) quickly deepens. The distance between 
the two peaks is in quantitative agreement with Ref. [1] . 
In addition, extensive comparisons with quantum Monte 
Carlo (QMC) have shown earlier [3,2] that our approach 
agrees quantitatively with QMC, and contains the same 
finite-size effects. In particular, p {uj — 0) is smaller in 
smaller lattices. Hence, while at T = 1/32 the crite- 
rion [l]p (oj = 0) < 1 x 10~ 2 is satisfied for N c = 64 

1/2 

and short £ ~ N c , we still need to verify that this 
reflects the behavior of a large (but not infinite) corre- 
lation length in the thermodynamic limit. That is why 
we verified that p(u = 0) < 1 x 10~ 2 for N c = 128 2 
as well, p (uj = 0) in dynamical cluster approximation 
(DCA) has the opposite size dependence and satisfies 
p(uj = 0) < 1 x 10~ 2 for N c = 64. Note that for size 
128 2 , £ already reaches 40 lattice spacings at T — 1/22. 
All of the above results may be understood analytically 
from the above equations [2] by considering the limit- 
ing case where the characteristic frequency in the spin 
spectral weight x" p becomes smaller than temperature, 
(renormalized classical regime). The local moment sum 
rule prevents a finite-temperature mean-field transition 
by letting U sp , and hence (n^ni), exhibit a downturn at 
T* . Below that temperature, £ grows rapidly but it be- 
comes infinite only at T — 0. Similarly, the opening of 




FIG. 1. Density of states as a function of u for various 
lattice sizes, the largest /3 having the largest peak heights. 
Inset shows (nfnj.) with symbols from Ref. [1], except for x 
and lines that represent our calculation. 

the pseudogap with decreasing temperature can be traced 
[2], in al — 2, to the singular contribution of Xsp (?) to 
£ct (fc) when £ becomes larger than the single-particle 
thermal de Broglie wave length £ t h = vf/T. Indeed, in 
that limit, the single-particle spectral weight A(k.H,oj) 
at hot spots is given by -2E" [(w - S') 2 + S" 2 ] _1 with 
lu - £' = and E" (k H , 0) oc f- d /( th - Since £/&h grows 
exponentially in the d = 2 renormalized classical regime, 
A (k# , oj) can become exponentially small at uj = even 
without a MIT. In the analytical approach, [2] the down- 
turn in (n-|-nj.) and the opening of a deep pseudogap are 
both unambiguously driven by a rapidly growing £ in the 
SDW channel. The pseudogap itself is not needed to re- 
inforce the downturn in (n^ni). While the situation is 
more subtle than that of Slater, the peaks in Fig.l are 
precursors of the SDW insulator that appears at exactly 
T = by the Slater mechanism. The peak separation in 
frequency (the gap) is larger than T* because Kanamori 
screening strongly renormalizes T* down. Increasing N c 
in DCA effectively lowers the dimension towards d = 2, 
revealing the effect of £ > £th on S and p. 
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